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Abstract. We revisit the calcnlation of dark matter relic abnndances in scalar-tensor gravity 
nsing a generic form for the conpling between the scalar held and the 

metric, for which detailed Big Bang Nncleosynthesis constraints are available. We hnd that 
BBN constraints restrict the modihed expansion rate in these models to be almost degenerate 
with the standard expansion history at the time of dark matter deconpling. In this case the 
maximnm level of enhancement of the dark matter relic density was fonnd to be a factor of 
~ 3, several orders of magnitnde below that fonnd in previons investigations. 
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1 Introduction 

Combined observations of large scale structure, galaxy dynamics and the cosmic microwave 
background (among others) [1] now provide overwhelming evidence for the existence of non- 
baryonic dark matter with present density (68% C.L.) [2] 

= 0.1188 ±0.0010(1.1) 

where ^dm is the dark matter density as a fraction of the total mass-energy budget and 
h = 0.6774 ± 0.0046 is dehned by the present value of the Hubble constant Hq = 100 h 
km/s/Mpc. 

In spite of this evidence, dark matter particles have never been observed directly and a 
convincing description of their particle nature remains elusive. The data favour cold (non- 
relativistic) dark matter with WIMPs (Weakly Interacting Massive Particles) the most pop¬ 
ular theoretical candidate.^ 

A popular framework for the origin of dark matter is provided by the thermal relic sce¬ 
nario in which the dark matter particles are produced through thermal scatterings of back¬ 
ground particles in the cosmic bath (see e.g. [3]). At early times, when the temperature of the 
universe is high, frequent interactions keep the dark matter particles in equilibrium with the 

^Within this class lies the widely studied neutralino, the lightest supersymmetric particle in supersymmetric 
extensions of the Standard Model in which i?-parity is conserved. 
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background. As the universe expands and cools the dark matter interaction rate drops below 
the expansion rate, the particles fall out of equilibrium, creation and annihilation processes 
cease and the number density redshifts with expansion. The surviving ’relic’ particles consti¬ 
tute the dark matter density we observe today. This process is known as particle freeze-out 
and, for non-relativistic dark matter particles, typically occurs at Tf ~ m^/20 (where rriy^ is 
the dark matter particle mass). 

Whilst in equilibrium, the dark matter number density, n^, decays exponentially (i.e. 
~ ~ where T is the temperature of the universe) so that the present dark 

matter abundance depends sensitively on the timing of freeze-out: the longer the species re¬ 
mains in thermal contact with the background bath, the lower its density at freeze-out. In 
the standard cosmological model of cold dark matter (CDM) with a non-zero cosmological 
constant (A), denoted the ACDM model, particle freeze-out occurs during the radiation dom¬ 
inated era when the expansion rate H ~ T^/Mpi (where Mpi = 1.22 x 10^® GeV is the Planck 
mass). In this scenario, a dark matter species with a weak scale interaction cross section, 
a ~ Gpm^, freezes out with an abundance that matches the presently observed value (1.1). 
This is known as the ’WIMP miracle’ and strongly motivates thermal WIMP dark matter 
models. 

Despite the observational success of ACDM, current datasets leave the physics of the 
universe prior to Big Bang Nucleosynthesis (BBN) {t ~ 200 s) relatively unconstrained. If the 
universe experiences a non-standard expansion law at early times, and in particular during 
the era of dark matter decoupling, particle freeze-out may be accelerated (or delayed) and 
the relic abundance enhanced (or suppressed) [4-12] (see also [13]). 

In this article we investigate the relic abundance of dark matter in scalar-tensor theories 
of gravity [14-20] whose cosmological evolution deviates from the standard expansion history 
in the pre-BBN era. These theories are (in part) motivated by higher dimensional unification 
models, a la Kaluza-Klein [21, 22], where additional scalar fields arise through the compact- 
ification of the higher dimensions and couple to the metric with gravitational strength. As 
such, the gravitational interaction is mediated by both the metric and scalar fields so that 
scalar-tensor gravity models represent a departure from standard General Relativity (GR). 

Importantly, the new long range interaction introduced by the coupling, between 

the scalar field, (/?*, and the metric, is subject to stringent experimental bounds from 

fifth force searches and solar system tests of gravity [23]. To evade detection in high density 
environments non-linear effects act to shield the scalar field through one of several screening 
mechanisms such as the Vainshtein or chameleon mechanisms [24].^ In the chameleon case 
for example, the mass of the scalar field is background dependent so that in regions of high 
density (e.g. our solar system) the field mass is large and the interaction range is suppressed. 
Gonversely, in low density backgrounds (i.e. on cosmological scales) the field can be extremely 
light, allowing the scalar interaction to drive the accelerated expansion of the universe. 

In addition to the screening mechanisms that help shield the scalar field from astrophysi- 
cal observations, many scalar-tensor gravity models exhibit an inherent attraction mechanism 
towards General Relativity [26, 27]. Throughout its cosmological evolution, the scalar field is 
driven towards a state where the coupling A{(p^) remains constant so that the scalar-tensor 
theory is indistinguishable from GR (see section 4). This allows scalar-tensor gravity models 

■^Matter fields ^'mat couple directly to the so-called Jordan frame metric where an 

asterisk is used to denote Einstein frame quantities, the frame in which the gravitational field equations take 
their standard form (see later). 

^For a recent paper on laboratory searches for the chameleon (scalar) field see [25]. 
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to deviate from the standard cosmological scenario at early times whilst relaxing towards 
General Relativity prior to the onset of Big Bang Nncleosynthesis. Hence, these models conld 
potentially distnrb the timing of dark matter deconpling and, in tnrn, the predicted dark 
matter density. 

The effects of the modified early time expansion rate in the scalar-tensor scenario on 
the dark matter relic abnndance were first stndied by Catena et al in [28]. In their paper the 
anthors determined the evolution of the scalar field for the conpling = 1 + and 

fonnd that, for T greater than some transition temperatnre the early time expansion 
rate was enhanced by a factor of 

^(T) = ~ 2.19 X 10^^ (where Tq = 2.35 x 10“^^ GeV) (1.2) 

^GR / 


before rapidly dropping to 1 at T = (as a resnlt of the attraction mechanism jnst men¬ 
tioned). Interestingly, the anthors discovered that the rapid relaxation of the scalar-tensor 
expansion rate, i/sT) towards the standard expansion rate, at Isd to a phase of 

reannihilation: after the initial particle deconpling, the dark matter species experienced a 
snbsidiary period of annihilation as the expansion rate of the nniverse dropped below the 
interaction rate. Despite this secondary annihilation phase, they showed that the relic abnn¬ 
dance of dark matter in scalar-tensor gravity models can still be enhanced by np to three 
orders of magnitnde (depending on the WIMP mass rn^). 

The approximate form (1.2) of the ratio has since been adopted in several 

snbseqnent investigations, inclnding those by Gelmini et al [9], Rehagen and Gelmini [29] and, 
very recently, Wang et al [30] who stndied the relic abnndance of asymmetric dark matter 
species [9, 30] and sterile nentrinos [29], and natnrally obtained similarly large enhancement 
factors for the dark matter relic abnndance. Snbseqnently Gatena et al [5] considered the 
conpling H((/?*) = 1 -|- as part of a more general stndy of dark matter relic abnndances 
in non-standard cosmological scenarios and once again fonnd that the relic abnndance in 
scalar-tensor cosmology is enhanced by several orders of magnitnde. 

We shonld emphasize, however, that althongh the different conplings nsed by the Gatena 
gronp allow for significant enhancements of the early time expansion rate, detailed stndies of 
the BBN implications of scalar-tensor theories with these conplings are lacking and only the 
simple constraint 


A((y7^,)|T=iMeV 

A((^*)|t=To 


(1.3) 


was imposed.^ However, detailed BBN stndies are available (see for instance [31, 32]) for the 
more popnlar choice of the scalar conpling 




(1.4) 


Thns it is important to determine if the relic abnndance of dark matter can be similarly 
enhanced for these more widely investigated models. 

The qnadratic conpling (1.4) is the simplest generalization of the original Jordan— 
Fierz—Brans—Dicke model (see later). Besides being the prototypical conpling considered in 
nnmerons investigations of scalar-tensor gravity theories, inclnding stndies of slowly rotating 

^This constraint on the coupling at the time of Big Bang Nucleosynthesis (T = 1 MeV) was derived from 
a constraint on the number of additional neutrino species AN < 1. 
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anisotropic neutron stars [33] and pulsars in binary systems [34-36], the quadratic coupling 
was recently the subject of the rigorous investigation by Coc et al [37] in which BBN calcu¬ 
lations were used to place stringent constraints on the various coupling parameters. These 
constraints are up to several orders of magnitude stronger than the solar system bounds 
(e.g. [38]) and, as we will show in the following, severely limit the possible deviations from 
the standard cosmological expansion history during the era of dark matter decoupling. 

Finally, we mention the 2008 paper by Catena et al [39] where the authors introduced 
a new hidden matter sector that experienced a different coupling than the ordinary visible 
sector, thus allowing for slower pre-BBN expansion rates. This non-universal scalar-tensor 
theory was found to predict dark matter relic abundances that were reduced by up to a factor 
of ~ 0.05. 

In this article we revisit the calculation of dark matter relic abundances in scalar-tensor 
gravity from first principles. We assume that the coupling between the scalar field and matter 
(including dark matter) is universal; the case of non-universal scalar-tensor theories, in which 
the coupling between the scalar field and the visible and dark matter sectors is distinct, 
will be briefly considered towards the end of the paper. We explicitly determine the modified 
expansion rate by solving the equation of motion for the scalar field directly and use this result 
to numerically solve the Boltzmann rate equation governing the evolution of the dark matter 
number density. We find, for the coupling A{Lp^) = that the efficiency of the attraction 

mechanism towards GR combined with the strict BBN bounds on the input parameters only 
allow for modest enhancements of the dark matter relic abundance, particularly compared to 
those reported in Catena et al [5, 28]. 

To begin, in section 2 we discuss the formulation of scalar-tensor theories in different 
conformal reference frames and comment on their physical interpretation. Then, after deriv¬ 
ing the equations that govern the cosmological evolution in section 3, we explore in detail the 
dynamics of the coupled scalar field and the attraction mechanism towards General Relativity 
(section 4). We then investigate which regions of parameter space lead to the largest devia¬ 
tions from the standard expansion history (section 5) and, most importantly, which regions 
satisfy the bounds from Big Bang Nucleosynthesis and other astrophysical and dynamical 
constraints (section 6). In sections 7 and 8 we calculate the dark matter relic abundance for 
symmetric and asymmetric dark matter models respectively, and, for the first time, determine 
the annihilation cross section required to produce the observed dark matter density. Then, in 
section 9, we consider non-universal scalar-tensor theories and discuss the issues associated 
with determining the dark matter relic abundance when the scalar interaction with the stan¬ 
dard matter (or visible) particle sector and dark matter particle sector is distinct. Finally 
in section 10 we summarize our results and comment on the potential for relic abundance 
calculations to discriminate between the scalar-tensor and standard cosmological scenarios. 

2 Scalar-tensor gravity: Jordan and Einstein frames 

Scalar-tensor gravity theories are often formulated in one of two conformal frames of reference, 
namely, the Jordan and Einstein frames.® The advantage of using these different conformal 

® Conformal reference frames are those connected by a conformal transformation, i.e. a local rescaling of 
the metric: 

( 2 . 1 ) 

Here, Q,{x) is an arbitrary function of the spacetime coordinates x^, and we use the notation to preserve 
the sign of the line element ds^ = OJdsi. 
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frames is that the scalar coupling enters through either the gravitational sector (Jordan frame) 
or the matter sector (Einstein frame), leaving the other sector unaffected. For example, the 
general action integral for a scalar-tensor theory with a non-minimally coupled scalar field 
formulated in the Jordan conformal frame, is given by 

5tot = j - 2Ui^)] + (2.2) 

where F(ip), Z{ip) and U{ip) are arbitrary functions of the field, G* is the bare gravitational 
constant (i.e. the gravitational constant in the absence of the scalar interaction), and R = 
g^’^R^v is the Ricci scalar which has been constructed from the Jordan frame metric g^y. The 
original Brans-Dicke model corresponds to F[(p) = Lp, Z{(p) = a;/(/? (where w is a constant) 
and U{ip) = 0. 

In this frame matter fields, ikmat, are coupled directly to the metric, gfj,y, so that the 
Weak Equivalence Principle is preserved by construction. This means that observables such 
as mass, length and time take their standard interpretation in the Jordan frame, making 
it the most convenient for particle physics considerations. However, since in this frame the 
scalar field couples to the gravitational sector, gravitational couplings pick up an additional 
(p dependence and the field equations take the cumbersome form: 


F{ip)G^y - - g^yU) F{ip) = SttG* {T^,y + , (2.3) 

where the covariant derivative V^, □ = and Einstein’s tensor, G/j^y = R^y — ^g^yR, are 

defined using the Jordan frame metric g^y, and the matter and scalar field energy-momentum 
tensors are given respectively by 


T — 

^ fiiy — 


2 


(2.4) 


and 


SnG^T^Py = Z{ip)dfj,ipdy<p - g^y 


]^g°‘dz{ip)da^dpip + U{ip) 


(2.5) 


To simplify the gravitational field equations we could formulate the scalar-tensor theory 
in the ’Einstein’ frame, which is related to the Jordan frame by the conformal transformation 


g^u = F{‘f)gi^u- ( 2 - 6 ) 

Here (and in the following) we use an asterisk to distinguish Einstein frame quantities from 
their Jordan frame counterparts. Applying the transformation (2.6), the action (2.2) becomes 


•S'tot = J d'^x [R* - d^tp^dyip^ - FV{(p^)] + 5mat[-4^((^*)5))v; ^mat] (2.7) 

where we have made the following redefinitions 


[ dip^'' 


d In F(ip) 

V dip , 

I 4 

dip 


Z{^) 

2F{ipy 


A{ipy=F-^/\ip), 


2v{ipy 


uM. 

F^ipY 
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The gravitational sector is now of the same form as a minimally coupled quintessence 
model. Moreover, the held equations take the simplihed form 


G;, = 87rG, (r;, + r;;) 

where the matter energy momentum tensor is now given by 

rp* _ 2 hS'mat 

and the scalar held energy momentum tensor is 


AttG^T^^* = - g 


fliy 


,g'^^daP*dpp^ + v{ip^ 


(2.9) 


( 2 . 10 ) 


( 2 . 11 ) 


However, in the Einstein frame, the scalar coupling enters through the matter action S'mat = 
Sma,t[-^^{'P*)9iiu', 'kjnat] SO that the matter helds Tmat couple to the (^*-dependent metric 
This indicates that particle physics quantities (e.g. mass, length, time, energies, 
cross sections) measured in this frame are spacetime dependent. 

We can characterize the departure of scalar-tensor theories from General Relativity by 
introducing the parameter 

d\nA{p^) 

«(</?*) =- -1 -• (2-12) 

cLiPif^ 

Large values of a{LpGj correspond to large variations in the coupling A({p^) whereas in the 
limit a(</?*) —>■ 0, corresponding to H((/?*) = const., the Einstein and Jordan frames coincide 
and the scalar-tensor theory reduces to standard General Relativity. 

In the following our strategy will be to determine the cosmological evolution in the 
Einstein frame, where the cosmological equations take their simplest form, and then transform 
our results over to the Jordan frame, which is where we solve the Boltzmann equation for the 
dark matter number density. 


3 Cosmological equations 

If we introduce the Einstein frame line element 

dsl = g^iydxi^dx'^ = -dtl + al{t)'y*jdxldxi, (3.1) 

and assume that the matter helds are a perfect huid described by the usual energy-momentum 
tensor 

T;, = {p*+p*)u;u:+p,g;,, ( 3 . 2 ) 

where p* and p* are the Einstein frame energy density and pressure of the huid respectively, 
then the held equations (2.9) give 

= SttG^p^ T 2V(p*), (3.3) 

3 — = — AttG^ (p* + p*) — 2(p^ -|- 2R (p*)- (3-4) 

Here an overdot denotes differentiation with respect to G and we have introduced the Einstein 
frame Hubble factor JL* = dlna^/dt^. Since the Friedmann equations (3.3) and (3.4) depend 
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on the scalar field, we also need the eqnation of motion for to close the system and 
determine the evolntion of the scale factor Variation of the action (2.7) with respect 

to gives® 

dV 

+ -— = —A'KG^a{(p^) (/?* — 3p*). (3.5) 

Cjipif. 

Additionally, the Einstein frame energy density and pressnre no longer satisfy the standard 
continnity eqnation. Instead there is an additional sonrce term dne to the scalar field inter¬ 
action: 

+ 3H^ {p^i + p^i) = a{(p^) {p^i - 3p*i), (3.6) 

where i labels the varions flnid components (e.g. radiation, baryonic matter, dark matter). 
Instead of solving this eqnation directly we can transform to the Jordan frame in which 

dpi 


dt 


+ 3H {pi + Pi) — 0, 


(3.7) 


where the Jordan frame expansion rate, H = dlna/dt is related to the Einstein frame expan¬ 
sion rate via 

H = + a{pG}p*\- (3.8) 

Note that (3.7) can be derived from (3.6) by transforming to the Jordan frame coordinates 

a = A(v9*)a* , dt = A{ip^,)dU, (3.9) 

and snbstitnting in the relationships for the Einstein frame energy density, p*, and pressnre, 

p* = A'^{(p^)p, p* = A^((p*)p. (3.10) 

Assnming the pressnre and energy density of the i-th flnid component are related by 
Pi = WiPi, where Wi is the eqnation of state parameter, we get 


Pi{t) = Pi (to) exp 

Farther, if we assnme that Wi is constant. 


Pi = Pi 


—3 / + Wi)dt 


'to 


-3(1-1-tOi) 


(3.11) 


(3.12) 


where a (snb)snperscript ’0’ denotes a quantity evaluated at the present epoch, i.e. p® = Pi(to)- 
Then, nsing (3.10), we finally have 


- P: 


0 


A{(p^ 

A{ip 




4-3(l+uii) / \ -3(l+uii) 

0*0 


(3.13) 


where p®j = A^((p*o)/Oi • Notice from the relationship (3.10) between the Einstein and Jordan 
frame energy densities and pressnres that the eqnation of state parameter is a frame invariant 
qnantity: 

(3.14) 


Pi P*i 

Wi = — = -. 

pi P^i 


®Also needed is the relation 


Stp* 
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4 Scalar field dynamics 


To determine the expansion rate of the universe in scalar-tensor gravity models we must 
solve the coupled system of equations (3.3)-(3.5). Fortunately, as pointed out by [26, 27], the 
equation of motion for the scalar field can be decoupled from the system by transforming the 
evolution parameter from the Einstein frame time to N = ln(a*/a*o)- In this case variables 
of the type Q transform as Q = Q '(where a ' denotes differentiation with respect to N)^ 
so that (3.3) becomes 

(3 - -h 2y((^*). (4.1) 

Since the potential term only affects the cosmological evolution at late times and does not play 
a significant role during the early universe period in which we are interested, i.e. E(y 5 ^,) <C 
~ p*, we set E((p*) = 0.^ Then, substituting (4.1) into (3.5) we get 

o ^ ,2 ^* + (1 “ - 3n^) 

3 - (p'Z 

where w is the equation of state parameter of the total cosmic fluid: 

^ _ P*tot _ Ptot 

P*tot Ptot 

Following the analogy given in [27], the field 99 * can be thought of as a particle-like dynamical 
variable with a velocity dependent mass, m((p*) = 2/(3 — <p(?). In this instance the parti¬ 
cle (^ 5 ^,) experiences simple damping and rolls down the potential oc ln4((p*) towards the 
minimum, provided such a point exists. Hence, the late-time evolution of the field is reason¬ 
ably straightforward given that the equation of state parameter during the matter dominated 
epoch is u) ~ 0 so that the forcing term on the right hand side of (4.2) is simply given by 
—a((p*). Therefore, if the function a((p*) possesses a zero with a positive slope, the field will 
be dynamically attracted towards the point a = 0, which is precisely the GR limit. For the 
choice of coupling 4((p*) = considered here, this condition implies that /3 > 0 . 

At much earlier times, deep the in radiation era {w ~ 1/3), the forcing term in (4.2) 
is mostly ineffective and any initial velocity possessed by the field is rapidly damped away. 
This allows us to take as initial conditions 

(p*in = const, and = 0, (4.4) 

where (p*in and are the values of the field and its A^-derivative at some initial point N[^. 

Although l — 3w is ~ 0 throughout most of the early universe, there is an important effect 
that arises when each of the particle species in the cosmic fluid becomes non-relativistic [27]. 

^The choice F(</5*) = 0 does not give a viable late time cosmology in that it cannot account for the recent 
accelerated expansion of our universe. Coc et al [37] have studied the effects of including a cosmological 
constant (in either the Jordan or Einstein frames) tuned to the observed value of the cosmological constant 
density today and find that only the late time dynamics of the scalar field is affected. Moreover, their analysis 
shows that the constraints on the parameters (/3, |ao|) (see section 6) are only moderately affected by the 
inclusion of a non-zero potential term, and that adding a cosmological constant in the Jordan frame actually 
constricts the allowed region in the (/3, |ao|) plane. Alternatively, Mota and Winther [40] consider general 
chameleon models (of which the quadratic coupling considered here is a special case) and study the conditions 
under which the chameleon potential F(pt) possesses an attractor solution where the chameleon follows the 
minimum of the effective potential F(<p*) -I- pmA(p*). The authors found that the field will typically settle 
at this minimum before the onset of BBN provided ';$> H. However, a later study by Erickcek et al [41] 
concluded large changes to occur near the minimum, casting doubt on chameleon models in general. 


(4.2) 

(4.3) 
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As the temperature of the universe drops below the rest mass of each of the particle types 
they provide a non-zero contribution to the quantity 1 — 3w. This activates the forcing term 
in (4.2) and displaces, or ’kicks’ the held along the potential oc \nA{(p^). In this way the 
attraction mechanism towards GR is initiated during the very early moments of the universe 
(T ~ mt ~ O(IO^) GeV, where mt is the mass of the top quark), prior to the onset of Big 
Bang Nucleosynthesis. 

To explore this effect in more detail we hrst write the quantity 1 — 3w as 

1 o Ptot 3ptot 
1 —6W = - - 

ptot 

Ptot y ^ 


- 3pa] + Pra 


(4.5) 


where we have separated out the contribution from the relativistic particle species A and the 
non-relativistic particle species m (for which pm = 0). Assuming that the total energy density 
is dominated by relativistic matter and radiation during the crossing of each mass threshold, 
then Ptot — Pr = 7r^p*p(T)T'^/30 where p*p(r) is the total number of relativistic degrees of 
freedom. The kick function 


S(T) = 

A 


PA - 3pA 

ptot 


(4.6) 


takes the form 


s(r) = 

A 


9A 2 f 

9*p{T) 



± 1 


(4.7) 


where = thaIT and -l-(—) corresponds to fermions(bosons), since the energy density pA 
and pressure pA of each of the particles of type A are given by (see e.g. [3]) 


Pa{T) 

Pa{T) 



9A f 
6vr2 


exp (E/T) ± 1 
exp {E/T) ± 1 


E'^ dE, 
dE. 


(4.8) 

(4.9) 


Here qa is the number of internal (spin) degrees of freedom. It is important to note that the 
variable T in the expression for S(T) is the Jordan frame temperature, which we can relate 
to the parameter N through 


T[p,,N] = To 


A((p*o) 


9*siTo) 


A{p^) ig*s{T) 


1/3 




(4.10) 


Here Tq = 2.35 x 10“^^ GeV is the current temperature of the universe and (p*o = T*{To)- 
Using the Standard Model particle spectrum given in table 1, we numerically evalu¬ 
ate (4.7). The results for the contribution from each of the individual particle species (dashed 
curves) as well as their combined sum (solid black curve) are shown in hgure 1. The quantity 
S(r) has been evaluated in several other publications (e.g. [37, 42]) and in most cases the 
value of p*p(T) is assumed constant during the crossing of each particle threshold. How¬ 
ever, since p*p(T) actually decreases during this interval, these calculations underestimate 
the height of S(T), particularly for the hnal ’kick’. In our work we have maintained the 
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0.1 



10 '® 10 '“ 10 ’^ 10 '^ 10 "' 10 ° 10 ' 10 ^ 10 ^ 

T [GeV] 


Figure 1. Evolution of the quantity E(T) = 1 — 3w (solid black curve), defined in (4.7), over the 
radiation era. The contributions from each particle species listed in table 1 are also shown (coloured 
dashed curves). 



Particle 

Mass (GeV) 

gi 


t 

173.07 

12 


b 

4.18 

12 



1.78 

4 


c 

1.28 

12 

Fermions 


0.106 

4 

s 

0.095 

12 


d 

4.8 X 10"3 

12 


u 

2.3 X 10-3 

12 



5.11 X 10“^ 

4 


V 

- 

6 


H 

125.00 

1 


Z 

91.19 

3 



80.39 

6 

Bosons 


0.140 

1 



0.135 

2 


7 

- 

2 


9 

- 

16 


Table 1. Properties of the Standard Model particle spectrum including mass (GeV) and number of 
spin degrees of freedom gi used in the evaluation of E(T'). Composite particles are marked with an 
asterisk. 


temperature dependence of g^p{T) and our results agree with those contained in the recent 
paper by Erickcek et al [43]. 

After the final kick, the forcing term remains inactive until the universe transitions 
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from the radiation dominated era to the matter dominated era in which 


1 - 3w 


Pm 


Pm + Pr 1 + T/Tf 


eq 


(4.11) 


where T^q ~ 0(10“®) GeV is the temperature of matter-radiation equality, i.e. p^iT^q) = 
Pm(Teq)- Combining the results shown in hgure 1 with the late time behaviour described by 
equation (4.11), we obtain the evolution of 1 — 3w from the radiation era up to the present 
time. This is shown in hgure 2 where we observe the general features that, in the radiation 
era, T > Tgq, 1 — 3w is ~ 0 except at the location of each of the kicks; then, as T approaches 
Teq, 1 — 3w smoothly rises to 1/2 and, in the limit T <C Tgq, 1 — 3w is ~ 1. 



Figure 2. Evolution of the quantity E(T) = 1 — 3w defined in (4.7) from the radiation era through 
to the matter era. 

The cosmological evolution of the scalar field can now be determined from (4.2). 
Substituting the computed values of 1 — 3w and our particular choice of coupling A{(p^) = 
into (4.2), we numerically integrate the equation of motion for several sample values of 
j3 and (^*in (see figure 3). 

In each case we see that the field is attracted towards y?* = 0 which, for the quadratic 
coupling considered, corresponds to the GR limit. Importantly we notice that the attraction 
mechanism is initiated well before the radiation-matter transition at T 10 ® GeV due to 
the non-relativistic kicks mentioned above. In fact, for the /? = 1 cases (blue and red curves) 
we can discern the four distinct kicks corresponding to the four peaks in figure 1. For the 
(/3, <y9*in) = (5,1) case (yellow curve), the attraction mechanism is so efficient that the field 
approaches the GR limit at much earlier times, prior to BBN. 

5 Modified expansion rate 

Having calculated the cosmological evolution of we can now determine the modified ex¬ 
pansion rate in the scalar-tensor scenario. For the purpose of calculating the dark matter relic 
abundance, we are particularly interested in the Jordan frame expansion rate H = dlna/dt 
since the dark matter particles couple directly to the Jordan frame metric and in this 
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Figure 3. Evolution of yi* for different values of the input parameters /3 and 


frame particle masses, number densities, etc. take their standard form. Therefore the ex¬ 
pansion rate that governs the timing of particle decoupling and in turn the dark matter relic 
abundance is determined using (3.8): 

H = + (5.1) 

where, using (4.1) and (3.10), we can write the Einstein frame expansion rate, ff*, as 

( 5 - 2 ) 

3 - 9 ?'/ 

To compare the Jordan frame expansion rate, H, to the expansion rate in the standard 
cosmological scenario, 

rr2 _ SttG 

we note that the gravitational couplings used in each case are related via [18] 

G = [l -|- ) (5-4) 

where G is the gravitational coupling given in the standard General Relativity scenario and G* 
is the bare gravitational coupling used here. Combining this expression with the connection 
between the Jordan and Einstein frame Hubble factors (3.8) we hnally have 

Hst _ A(ip^) l + a((p^)(p'^ 1 

Hqr sjl- Y'T+~a2(^^ 

where we now denote the Jordan frame expansion rate H, that will be used as input into the 
Boltzmann rate equation for the dark matter number density, by i^sT- 

To estimate the possible enhancement of the dark matter relic abundance in the scalar- 
tensor gravity scenario, we hrst determine the magnitude of the ratio ^ = Hs'j'/Hqh around 
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Figure 4. Magnitude of the ratio of the scalar-tensor and standard expansion rates, ^ = Hst/H an, 
evaluated at a temperature T = 5 GeV as a function of the input parameters j3 and which satisfy 
the dynamical constraints. 


the time of dark matter decoupling. Hence, in figure 4, we plot the magnitude of ^ evaluated 
at Tf ~ my^j20 ~ 5 GeV, which is the typical freeze-out time for a 100 GeV WIMP. We have 
chosen values of /3 > 0 to ensure that the model is dynamically attracted towards GR (see 
previous section) and omitted data points for larger /3 and that violate the dynamical 
constraints, leading to unphysical results (see next section). 

This figure shows that the scalar-tensor expansion rate at the time of dark matter de¬ 
coupling can be more than an order of magnitude larger than the standard expansion rate. In 
particular, the largest enhancements are seen for small values of (5 and large specifically, 
for /3 < 2 and > 2. This parameter range leads to large initial values of ^in ~ ^(<P*in) 
(because of the large (^*in) whilst exhibiting a less efficient attraction towards GR (small (5) so 
that the expansion rate at the time of freeze-out still deviates significantly from the standard 
one. In general, although increasing both /3 and increases the initial value ^in, because 
the displacement of the field towards y?* = 0 due to each of the non-relativistic ’kicks’ also 
increases with increasing /3 and the overall effect on ^(T/) (and in turn on 
can be difficult to predict. 

Although the early time expansion rate can certainly be much larger in the scalar-tensor 
scenario for certain regions of parameter space, we must be careful to ensure that such regions 
satisfy both the astrophysical and dynamical constraints placed on the scalar field and its 
evolution. In the next section we investigate in more detail these constraints and determine 
those points in parameter space which are viable. 

6 Constraints 

Solar system tests of gravity, including the perihelion shift of Mercury and Lunar Laser 
Ranging experiments, place strict constraints on deviations from General Relativity (see for 
example [23]). Most relevant for our purposes are the measurements of the Shapiro time delay 
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performed by the Cassini spacecraft [38] which indicate that the present value of the scalar- 
tensor deviation parameter Oq = q;^((^*o) < 10“^. Hence, if the gravitational interaction is 
truly described by a scalar-tensor theory, it must be extremely close to General Relativity in 
our local neighbourhood.® 

Although this constraint only applies to the present value |ao|) we can relate it to the 
input parameters (5 and by integrating the evolution equation up to the present epoch 
and calculating the predicted value of |ao| directly. We can then determine which values of 
j3 and satisfy the Cassini bound Oq < 10“^. The results are shown in figure 5 where we 
have indicated those points that violate the Cassini bound with a red cross and those that 
are acceptable with a green circle. We have also indicated with a blue cross those points that 
violate the various dynamical constraints on the evolution of the scalar field (see later). 



Figure 5. Exclusion plot following a scan of the input parameter space (/?, y>*in) and applying the 
Cassini bound Oq < 10“® [38]. The points excluded due to the Cassini (dynamical) constraints have 
been marked with a red (blue) cross and the allowed points are indicated by the green circles. 


Interestingly, for the parameter range considered, the Cassini bound only excludes values 
of /3 < 0.4. Comparing this with the results in figure 4 for ^(Tf), we see that the regions of 
parameter space that lead to the fastest expansion rates at the time of dark matter decoupling 
(0.4 < /3 < 2,(^*in > 2) still satisfy the constraint Og < 10“^. Therefore, the Cassini bound 
alone (or, more generally, solar system tests of gravity) does not preclude large deviations 
from the standard cosmological history at early times. This is not so surprising given that 
solar system tests of gravity rely on data taken at late times — long after the attraction 
mechanism towards GR has been initiated. 

To properly constrain the input parameters and the behaviour of the field in the early 
universe, we must also take into account the results presented in Coc et al [37] where they con¬ 
sidered the Big Bang Nucleosynthesis implications of scalar-tensor theories with a quadratic 
coupling, i.e. In this work the authors perform a full numerical integration of 

®Additionally, for the quadratic scalar coupling A(ip*) = considered here, the decay of the orbital 

period of pulsars in asymmetric binaries implies that the coupling parameter j3 > —4.5 [36]. 
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the scalar field evolution equation and calculate the resulting light element abundances using 
an up-to-date BBN code. The results of these calculations are then compared to observed 
light element abundances to constrain the various scalar-tensor model parameters. 

Since these constraints are also given in terms of the present value |ao|) we must again 
integrate the evolution equation up to the present epoch for the different values of /3 and 
and calculate |ao| directly. As an example, in figure 6, we compare our results for |ao| 
as a function of /? for different starting values of (solid curves) with both the Cassini 
(dot-dashed purple curve) and the Coc et al BBN bound (dot-dashed black curve). 



Figure 6. Magnitude of jaol as a function of the coupling parameter /3 for different initial values of 
the field (solid curves). The results are compared with the constraints derived using the solar system 
tests of gravity (dot-dashed purple) and the BBN constraints given in figure 19 of Coc et al [37] 
(dot-dashed black). 

Ignoring any computational differences, (keep in mind that our computation of S(T) 
differs from figure 5 of Coc et al [37] because we maintained the temperature dependence of 
g^p{T)) we notice that the blue, red and yellow curves for jaol do not drop below the BBN 
constraints (dot-dashed black curve) until after a particular oscillation. For example, for 
<y9*in = 0.75 (red curve), the calculated values of |ao| do not drop below the BBN constraints 
until after the fourth oscillation at ,0 ~ 1.65. Similarly, for = 1 (yellow curve), |ao| does 
not drop below the BBN bound until after the fifth oscillation at /3 = 2.35. Therefore, once 
the BBN bounds become more stringent than the solar system bound, i.e. /? > 0.4, the values 
of (5 for which a particular becomes acceptable are discrete. This is shown in figure 7 
where we reproduce figure 5, this time applying the Coc et al BBN bound in addition to the 
Cassini bound. For /3 < 0.4, the boundary separating the allowed and rejected regions follows 
a smooth curve, whilst for /3 > 0.4 the boundary increases in discrete steps. 

Most importantly, a comparison of figures 5 and 7 shows that applying the Coc et al 
BBN constraints excludes most of the parameter space corresponding to large early time 
expansion rates (see figure 4). Also, because of the non-injective relationship between the 
input parameters and the present value of jaol (be. different initial conditions can lead to the 
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Figure 7. Scan of the input parameter space for the coupling A{ip^) = The points excluded 

due to astrophysical (dynamical) constraints have been marked with a red (blue) cross and the allowed 
points are indicated by the green circles. 


same |ao|)) the constraints in figure 19 of [37] can only be treated as a strict upper bound. 
This means that values of |ao| lower than this upper bound are not necessarily acceptable so 
that the allowed region in the (/3, ip^n) parameter space is most likely smaller than what is 
shown in hgure 7. 

In addition to the astrophysical constraints, the dynamical evolution of the scalar held 
must also satisfy the condition 


1 + > 0 . 


( 6 . 1 ) 


This ensures that the Jordan frame Hubble factor, derived from the Einstein frame Hubble 
factor through 

H = [l + , (6.2) 

is positive dehnite. Moreover, if this condition is violated, the Jordan frame temperature T 
becomes a multi-valued function of the evolution parameter N and the conformal transfor¬ 
mation between the Jordan and Einstein frames breaks down.® The points in the (/3, 
plane that violate this dynamical constraint are indicated with blue crosses. 

We can crudely estimate the general shape of the dynamical exclusion zone by approxi¬ 
mating the maximum velocity reached by the held due to each ’kick’ as ~ ck( 7’*)- The 


®If we neglect the temperature dependence of gr,s(T), the connection between the Jordan frame temperature 
and the evolution parameter N (4.10) gives 


dT 

dN 


A{ip,o) _jv 

A{p.) ® 


[1 -b a(y>*)‘P*] • 


(6.3) 


If the quantity 1 -|- a((^*)(^(, changes sign, the relationship between T and N is not monotonic and T becomes 
a multi-valued function of N that cannot be inverted. 
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condition 1 + a{ip^)^p'^ > 0 then becomes 


a{^.f<V\ (6.4) 

where T) is some variable that depends on the magnitnde of the kick. Snbstitnting in the 
expression «((/?*) = (3^*^ we hnd that the valnes of in the region 

(6.5) 


are exclnded. 

Combining the astrophysical and dynamical constraints we see that a significant portion 
of the parameter space is exclnded. In particnlar, the region corresponding to the greatest 
deviations from the standard expansion history at the time of dark matter deconpling (i.e. 
CiTf)) is no longer allowed (compare hgnres 7 and 4). In fact the largest valne of C(^/) 
permitted by the varions constraints is only ~ 4. Hence, althongh large early-time expansion 
rates can be achieved within scalar-tensor cosmological models with the conpling H(y?*) = 
6 those models that satisfy the varions constraints mentioned above will not signihcantly 
deviate from the standard cosmological model at the time of dark matter deconpling. In 
tnrn, we expect that the relic abnndance of dark matter to be relatively nnaffected and that 
any enhancements with respect to the canonical resnlt are modest. In the next section we 
verify this assertion by nnmerically solving the Boltzmann rate eqnation governing the cosmic 
evolntion of the dark matter nnmber density to determine the present dark matter density. 


7 Symmetric dark matter 


The relic abnndance of a symmetric dark matter species x(= x), initially in eqnilibrinm with 
the backgronnd cosmic bath, is determined from the dark matter nnmber density which 
evolves according to the relativistic Boltzmann eqnation 

^ = -3Hn^ - {av) (nl - . (7.1) 

Here is the eqnilibrinm nnmber density and (av) is the thermally averaged annihilation 
cross section times relative velocity (loosely termed the "annihilation cross section"). We 
adopt the generic parameterization for the annihilation cross section [44, 45] 


(av) = a + 


6r 


(7.2) 


where the constant term, a, corresponds to s—wave scattering and the temperatnre dependent 
term, b, to p—wave scattering. 

The Boltzmann eqnation (7.1) can be rewritten in terms of x = m^/T and the comoving 
nnmber density Y = n^js, where s is the entropy density given by s = 2tt‘^ g^{T)T'^/ 45.^^ We 
then have 


dY 

dx 


s{av) 

xH 


C{x) (y2 


y2'| 

eqj 


(7.3) 


^°Here g, (T) actually refers to the number of entropic degrees of freedom . Since the number of relativistic 
and entropic degrees of freedom only differ when a particle crosses a mass threshold, we take g^p = g *3 = 
9* [46]. 
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where ~ Q.lAh{g^/g^,) g^. = 2 is the number of internal degrees of freedom of the 

dark matter species x and 



(7.4) 


is a temperature dependent factor related to the change in the number of degrees of freedom. 
The present dark matter density, is obtained from the asymptotic solution (x —>■ oo) 

of equation (7.3) 


= 2.75 X lO^m^yoo, 


(7.5) 


where Y^o = T(x —>■ oo) is the present comoving density. 

As discussed in section 2, each of the quantities in (7.3) such as the comoving number 
density, Y = n/s, and annihilation cross section, {crv), is dehned in the Jordan frame where 
particle physics quantities take their standard interpretation. Hence H is the Jordan frame 
expansion rate determined using (3.8).^^ 

Taking the results from our scan over the input parameter space, we numerically integrate 
the y?* evolution equation (4.2) and determine the Jordan frame expansion rate H from (5.1). 
In doing so we have chosen values of j3 and that are (i) permitted by the astrophysical 
and dynamical constraints (according to hgure 7) and (ii) provide the greatest deviation from 
the standard expansion history (see hgure 4). 

Substituting this result into (7.3) we solve the Boltzmann equation to determine the 
present dark matter density. The results are shown in hgure 8 where we plot the ratio of the 
dark matter relic abundance in the scalar-tensor gravity model, to the corresponding 

value in the standard cosmological scenario, (i-e. the solution to (7.3) using Hq^), as 

a function of for different values of /3. In the left and right panels of hgure 8 we have 
calculated the relic abundance enhancement ratio for = 10 GeV and rriy. = 100 

GeV respectively with the solid (dashed) curves corresponding to s(p)—wave annihilation. 
Note that for the values of (3 chosen, namely {3 = (2.4,3.0,4.0), we have only calculated 
results up to = (2.5, 2.0,1.5) respectively because larger values of (^*in are excluded due 
to dynamical constraints (see hgure 7). 

In each panel we see that the ratio increases for small reaching a (3- 

dependent maximum value before falling as increases further. For the = 10 GeV and 
= 100 GeV cases the maximum enhancement for s—wave annihilation is only ~ 

1.5 and 2.6 respectively and is obtained for (/3,(/9in) = (2.4,1.8) in both cases. For p—wave 
scattering and = 100 GeV, this ratio increases to ~ 2.9. 

Because the lack of an algebraic expression for the Jordan frame expansion rate H 
(see (5.5)) inhibits our ability to derive an approximate analytical solution for the dark matter 
relic density, it is difficult to make any quantitative estimates for the ratio in terms 

of (3 and Although we expect that the relic density enhancement factor, is 

related to the magnitude of ^ = F7 st/77gr evaluated at the time of dark matter decoupling, 
i.e. ^{Tf), the situation is still not straightforward given that the relationship between C{Tf) 
and the input parameters /3 and is difficult to predict (see the discussion in section 5). 
We simply comment that the relic abundance curves displayed in hgure 8 follow the same 
prohle as the corresponding ^(T/) curves evaluated at hxed 

the universal coupling case there is a unique Jordan frame in which both the Standard Model particles 
and the dark matter particles couple directly to the metric This situation becomes more complicated 

in non-universal scalar-tensor theories where the coupling with the visible and dark sectors is distinct (see 
section 9). 
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Figure 8. Ratio of the relic abundance of symmetric dark matter in the scalar-tensor and standard 
cosmological scenarios for = 10 GeV (left panel) and rrix = 100 GeV (right panel) as a function 
of (p*in for different values of /3. The solid (dashed) curves correspond to s(p)—wave annihilation. 

Most importantly, for the quadratic coupling A(<y9*) = e 2 ^‘^* considered here, the max¬ 
imum enhancement achieved for the s—wave case is only ~ 2.6 which, for the equivalent 
WIMP mass, is much less than the factor of ~ 400 reported by Catena et al [28]. Similarly, 
the maximum enhancement for the p—wave case, ~ 2.9, is several orders of magnitude below 
the factor of ~ 1400 given in [28]. 

Moreover, we find no evidence for reannihilation. That is, during our scan over the 
input parameters we did not observe a secondary phase of dark matter annihilation following 
particle freeze-out. In the Catena et al 2004 paper [28], reannihilation was a consequence of 
the rapid relaxation of the scalar-tensor expansion rate towards the standard expansion rate, 
which they purport occurs after dark matter decoupling. However, we showed in section 4 
that the attraction mechanism towards General Relativity is initiated when the temperature 
of the universe first drops below the rest masses of Standard Model particles at T > 10^ 
GeV. Even taking the most optimistic estimate that Xf ~ 10 in the scalar-tensor scenario, 
a dark matter particle with rriy- = 50 GeV (the same value used in figure 7 of [28]) would 
freeze out at Tj = 5 GeV. Therefore, by the time the dark matter particles decouple from the 
thermal background the attraction mechanism is typically well underway and the scalar-tensor 
expansion rate is already relatively close to the standard result. 

Using the observed dark matter density (1.1) we can invert our calculation to determine 
the required annihilation cross section. The results are plotted in figures 9 and 10 for the s— 
and p—wave annihilation cases respectively and show similar behaviour to the relic density 
curves given in figure 8. 

As in [11, 12] we could compare our results with the Fermi-LAT data [47] which places 
an upper limit on the dark matter annihilation cross section. However, since, in this case, the 
enhancement factors are so small, the Fermi-LAT bounds would not be able to significantly 
constrain the allowed model parameters. 

8 Asymmetric dark matter 

For completeness we also determine the relic abundance of asymmetric dark matter particles in 
scalar-tensor gravity. Since the symmetric annihilation cross section is only enhanced slightly, 
we expect a similarly small enhancement for the asymmetric annihilation cross section. 
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Figure 9. Required s—wave annihilation cross section {av) for symmetric dark matter for = 10 
GeV (left panel) and = 100 GeV (right panel) as a function of for different values of /?. 




Figure 10. Same as figure 9 but for p—wave annihilation. 


Asymmetric dark matter models treat the dark matter particle x antiparticle x 
as distinct and with unequal number densities, similar to the asymmetry that exists in the 
baryonic sector [48-50]. These models typically assume [51] either a primordial asymmetry 
in one sector that is transferred to the other sector, or that both asymmetries are generated 
by the same physical process such as the decay of a heavy particle. Relating the asymmetry 
in the dark matter sector to that in the baryonic sector also explains the proximity of the 
dark and baryonic densities, Q-oM/^b ~ 5, suggesting the dark matter mass is in the range 
~ 5 — 15 GeV [52]. 

For distinct particle x ^-nd antiparticle Xi the Boltzmann equation (7.1) is generalized 
to the coupled system 

— (av) {n^n^ — , (8.1a) 

3Hn^ — (av) {riy^n^ — , (8.1b) 

where and are the equilibrium number densities of the x x components respec¬ 
tively. We assume that self annihilations are forbidden, and that only interactions of the type 


dn^ 

dt 

dn^ 

dt 
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Figure 11. Iso-abundance contours in the {{av),C) plane for asymmetric dark matter correspond¬ 
ing to the observed dark matter abundance = 0.1188 for = 10 GeV (left panel) and 

= 100 GeV (right panel). The blue, red and yellow curves correspond to the input parameters 
(/3, (/3*in) = (2.4, 2.0), (3.0,2.0) and (4.0,1.5) respectively and the black curves correspond to the stan¬ 
dard cosmology result. We have calculated the results for both s— (solid curves) and p— (dashed 
curves) wave annihilation. 


XX (where the X’s are Standard Model particles) can change the dark matter particle 

number. We can then write 

= C, (8.2) 

where C is a strictly positive constant that characterizes the asymmetry between the parti¬ 
cles and antiparticles. Here, we are not concerned with the mechanism that generates the 
asymmetry, only that one has been created well before particle freeze-out. 

In terms of the comoving densities and Y^, the Boltzmann equations (8.1) become 


dYx 

dx 

dY^ 

dx 


s{av) 

xH 

s{av) 

xH 


C{x) {Y^ - CY^ - P) , 
C(x) (F| + CF^-P), 


(8.3) 


where, since the dark matter particles and antiparticles are non-relativistic at decoupling, 


P = F^'iFl'i 

X, X 


0.145 


9* 


x^e 


(8.4) 


Solving the system (8.3) in the asymptotic limit, the total dark matter density, is 

the sum of the x and x components, 


= 2.75 X 10® (F~ + Y^) . (8.5) 

In figure 11 we plot the iso-abundance contours in the {{(Jv), C) plane corresponding to 
the observed dark matter density = 0.1188. The results have been calculated for both 

= 10 GeV (left panel) and = 100 GeV (right panel) for both the s— (solid curves) 
and p— (dashed curves) wave cases. 

Because the vertical asymptote is only slightly shifted from the standard cosmology 
result (black curve), we do not expect the detection signal from annihilation of asymmetric 
dark matter to differ appreciably in the scalar-tensor gravity scenario. However, because 
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the curves are still shifted, it is possible to have a stronger asymmetric detection signal in 
the scalar-tensor scenario compared with the symmetric signal in the standard scenario (see 
e.g. [9, 11, 12]) . 

The relic abundance of asymmetric dark matter in scalar-tensor gravity has very recently 
been studied by Gelmini et al [9] and Wang et al [30]. However these studies are based 
upon the enhancement factor (1.2) and also neglect the temperature dependence of g*{T). 
Unsurprisingly, our results are much less dramatic than those given by Gelmini et al [9] and 
Wang et al [30] who found that the required annihilation cross section can be increased by 
several orders of magnitude. 


9 Non-universal scalar-tensor theories 
9.1 Visible and Dark Jordan frames 

The formalism and results derived thus far only apply for the special class of scalar-tensor 
theories known as universal scalar-tensor theories in which each of the different matter fields 
experiences the same coupling H((/?*). Here we consider a more general arrangement where 
the coupling to each matter sector is distinct. Specifically, we will consider the case studied 
by Goc et al 2009 [53] where the coupling between the scalar field and the visible (or Standard 
Model) sector, denoted Av{(p*), is distinct from that between the scalar field and the dark 
sector, In this case the Einstein frame action integral (2.7) is generalized to 

Stot = j - 4V((^*)] 

+ Sv[-Avi^*)9% i + ^D[-A‘jj{ip*)g*^i ^; 'I'^)], (9.1) 

where Sy and Sd are the action integrals for the visible, Ty, and dark, T/j, matter fields 
respectively. The corresponding generalization of the set of cosmological equations appearing 
in section 3 to the non-universal case is given in Goc et al [53]. We give a brief summary in 
Appendix A. Here we discuss the implications for the derivation of the Boltzmann equation. 

For the general case Ay ^ Ajj it is not possible to transform to a conformal frame in 
which both the visible and dark sectors couple directly to the same metric, that is, there is no 
unique definition of the Jordan frame. Therefore we must distinguish between two separate 
Jordan frames: the Visible Jordan Frame (VJF) defined by the metric 

gl, = Al{^,)gl^ (9.2) 

and the Dark Jordan Frame (DJF) defined by 

9^u = Aj{^*)9lu = B{ip*)gl^,. (9.3) 

In the VJF and DJF respectively, the visible and dark sector matter fields do not experi¬ 
ence the scalar coupling and particle properties (e.g. mass, energy, cross sections) and their 
interactions take their standard form. 

^^This set-up is different from that considered by Catena et al in [39] in which they introduced a new hidden 
matter sector and included dark matter within the visible sector. 
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9.2 Boltzmann equation 


Before attempting to calculate the dark matter relic abundance, we should first pause to 
contemplate whether, in the non-universal coupling case, the assumptions made in deriving 
the Boltzmann rate equation (7.1) (see e.g. [3]) remain valid. 

Starting with the fundamental form of the Boltzmann equation, 

i[/(x^P^)] = C[/(x^p^)] (9.4) 

where f{x^,p^) is the dark matter distribution function, L is the relativistic Liouville operator 
and C is the collision operator, the first question to address is in which frame should we 
formulate the problem? Given that the subject of the calculation, the dark matter particles, 
couple directly to the metric < 7 ^, the obvious answer is the DJF. Then, proceeding as in [3] 
we can determine the zeroth moment of (9.4) by integrating both sides of the equation over 
momentum space. 

To evaluate the integral of the Liouville operator, we note that the conformal transforma¬ 
tion (9.3) conserves both the isotropy and homogeneity of the VJF metric so that the dark 
matter distribution function, which is defined in the DJF, reduces to f{x^,p^) = 

Therefore, we recover the standard expression 

(i? /r|/(E.‘)l| = + 3HDno (9.6) 

where ud = nDitn) and to are the dark matter number density and time as defined in the 
DJF. 

However, the evaluation of the collision term in (9.4) is not so straightforward. The 
change in the dark matter particle number is governed by reactions of the type 

XX ^ (9.7) 

where the y’s are the dark matter particles and the X’s are particles belonging to the visible 
sector. In the non-universal case the particles on either the left or right hand side of the 
reaction (9.7) will either be susceptible or immune to the scalar interaction depending on the 
choice of conformal frame of reference; there is no conformal frame in which the scalar coupling 
vanishes for both sides of the reaction. Therefore, if we decide to formulate the Boltzmann 
equation in the DJF we must account for the varying (^^-dependent masses, energies, cross 
sections, etc. of the visible particles. Of course, transforming to the VJF does not alleviate 
the problem because then dark matter particles would be -dependent. 

A proper analysis of the problem would require a thorough reexamination of the deriva¬ 
tion of the Boltzmann equation and in particular the evaluation of the collision integral. Since 
this is beyond the scope of the present study we leave it as a suggestion for future work. 

^^Observations of the isotropy and homogeneity of the universe are made with Standard Model particles 
and therefore apply to the VJF. 

^"^In general, a conformal transformation is a local rescaling of the metric, 

= n'^{x)gf,„, (9.5) 

that only preserves the isotropy of spacetime and not homogeneity. But, if the conformal factor is a function 
of time only, i.e. = Q{t), as in (9.3), then homogeneity is conserved also. 
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10 Conclusions 


The inherent attraction mechanism exhibited by many scalar-tensor gravity models towards 
General Relativity allows for deviations from the standard cosmological scenario in the early 
universe that may not show up in present observational data. In fact, it has been conjectured 
that relic abundance calculations may be one of the few probes capable of discriminating 
the predictions of scalar-tensor scenarios from standard General Relativity [13]. To find out, 
we determined the evolution of the coupled scalar field from first principles, allowing us to 
calculate the modified expansion rate in the scalar-tensor gravity scenario, which we then 
used to calculate the dark matter relic abundance. 

As a specific example we considered the prototypical quadratic coupling A((p^) = 
and found that the maximum enhancement for a = 100 GeV WIMP was ^^dm/^dm ~ 
Although this ratio would increase with increasing WIMP mass, the level of enhancement is 
still far below that found in previous relic abundance investigations in scalar-tensor cosmol¬ 
ogy [5, 28]. 

Although the expansion rate at the time of dark matter decoupling in the scalar-tensor 
scenario can be up to several orders of magnitude larger than the expansion rate in the 
standard scenario (see figure 4), these points were excluded by BBN constraints [37]; had these 
points been acceptable, we would have found much larger relic density enhancement factors, 
possibly in line with those reported in [28] and [5]. Given that detailed BBN constraints for 
the couplings A{(p^) = 1 + and A{i~p^) = l + htpl, considered in [28] and [5] respectively, 

were not applied in their investigations because a study equivalent to [37] was not available, we 
suggest that the allowed expansion rate in those models may be much smaller than previously 
reported, and that the level of enhancement of the dark matter relic abundance might actually 
be much closer to the levels found here. 

To complete our study we also investigated the relic abundance of asymmetric dark 
matter species in scalar-tensor gravity models and found that the required annihilation cross 
section was enhanced only slightly, in contrast to the several order of magnitude estimates 
given in Gelmini et al [9] and Wang et al [30] who both used the parameterization (1.2) given 
in Gatena et al [28]. 

Finally, since the attraction mechanism towards GR is typically initiated prior to dark 
matter freeze-out so that the effect on the dark matter relic abundance is only modest, the 
scalar-tensor scenario with a quadratic coupling to matter does not provide a significantly 
different picture from the standard cosmological scenario. Hence, we conclude that unless 
reasonably precise details about the nature of the dark matter particle and its interactions 
are known, relic abundance calculations are most likely unable to discriminate scalar-tensor 
gravity from General Relativity. 


A Non-universal scalar-tensor theories: Cosmological equations 

The cosmological equations for a non-universal scalar-tensor theory for which the coupling 
between the visible (R) and dark (D) sectors is distinct, i.e. Ay ^ A^, has been given 
previously in, e.g. Goc et al [53]. We briefly summarize the main results here. 

Following the prescription given in section 2, we can apply the conformal transformation 

9^u = ^viv*)9lu (A.l) 
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to remove the coupling from the visible sector matter action so that the fields 'I'y couple to 
the metric directly. The action integral (9.1) then becomes 

•S'tot = J yj-gv [F{lp)Rv - g'^"'Z{(p)d^ipdu^ - 2U{lp)] 

+ Sv[g^y ; ^'y] + SD[B‘^{<f*)gl,y ; ^'d]; (A.2) 

where the connection between (p and and the definitions of the functions F{(p), Z{ip) and 
U{(p), are given in (2.8) with ^((/J*) replaced by Ay{(p^,). Also, we have defined 

Although the visible sector action Sylg^^ ; 'I'y] is independent of the scalar field in this new 
frame, the term SD[B'^{(pif)g'^y]'^D\ still contains a scalar coupling. We must therefore dis¬ 
tinguish between two distinct Jordan frames: the Visible Jordan Frame (VJF) defined by the 
metric g^^, and the Dark Jordan Frame (DJF) defined by 

= ADi^*)gly = B{ip^)gl^. (A.4) 


Similar to (2.12), we can define the parameters that characterize the deviation from General 
Relativity as 




din Aj((^*) 


(A.5) 


where j = V,D labels the different matter sectors. 

In the non-universal case the cosmological equations (3.3) and (3.4) for the evolution of 
the Einstein frame scale factor a* remain unaltered. However, the matter fields of each sector 
obey separate continuity equations 


dp^j 

dU 


J- 3id* {p*j F p*j) 


(XjippA) {^P*j 


(A.6) 


This implies that the visible (dark) sector matter fields are no longer conserved in the DJF 
(VJF) so that, for example, in the VJF we have 


PVi = Pvi 


\avoJ 


PD = Pd 


B{(p^) 

B{p*o) 


4—3(1+io_d) 


\ayoJ 


(A.7) 


where i labels the different visible fluid components and pn is the dark matter energy density 
as measured in the VJF. 

Moreover, the equation of motion (3.5) for the scalar field is generalized to 


F SH^ip^ -|- 


dV 

dp^ 


—JttG* ^ Oj (p^j — 3p^j) . 
j=V,D 


(A.8) 


Finally, the connection between the Einstein frame expansion rate and the VJF and DJF 
expansion rates is given by 

Hj = [id* aj(v9*)(^*] . (A.9) 
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Coc et al 2009 [53] have studied the BBN constraints on models with quadratic couplings 

{j = V,D). (A.IO) 

In particular, they studied regions in the {I3v,/3d) parameter space for which there is late 

time attraction towards GR and then derived additional constraints on the parameter space 

from BBN and precision gravitational tests. 
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